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Introduction 

During this phase of the program the analytical model for the amplifica- 
tion of broadband jet noise has been further developed and some numerical 
solutions have been performed. These calculations have been published [13 
(attached as Appendix A). The calculations on the stability of compressible 
axlsymmetrlc jets, which were described In the first annual report on this 
grant » have been documented and submitted for publication [2] (attached as 
Appendix B). The facilities. Instrumentation and data processing techniques, 
required for the experimental part of this program are nearing completion 
and some preliminary measurements have been performed. It Is expected, as 
described below, that much of the flow data, which is required to assist 
In the modeling of the flow will be collected during the second part of this 
grant. A brief discussion of the analytical and experimental parts of the 
program is given below. 

Analysis 

The present status of the calculations Is described in detail in 
references [1] and [2] so details are not presented here. However, some 
general conclusions can be drawn. The numerical procedure that has been 
developed is definitely more flexible than previous formulations. These 
previous works have employed an integral analysis which requires assumptions 
to be made about the characteristic shapes of the mean velocity profile and 
the turbulent kinetic energy distribution, (see Chan [3], for example). The 
present work confirms that these assumptions are reasonable in the case of 
a single Jet in a stationary ambient medium. The power of the present 
method lies in its ability to examine interactions in flows with more 



complicated geometries. For example the present technique Is readily 
adaptable to the problem of a single jet In a co-f lowing stream with a 
significant boundary layer on the external surface of the nozzle. Such 
calculations were not possible with the integral formulation. During the 
second half of this year> while awaiting experimental results, work will 
concentrate on improving the calculation procedure. Specific attention will 
be paid to the initial conditions. The present calculations have assumed 
a mean velocity profile at the jet exit which Is already developed Into the 
similarity form of the annular mixing region. Calculations will be performed 
astmmlng &, turbulent boundary layer at the jet exit. The effect of initial 
boundary layer thickness on the pure tone excitation of the jet will then 
be examined. The model proposed for the "wave-induced stresses" In the 
present calculations involves a scalar eddy viscosity. Until the experi- 
mental results are available no further development of this model will be 
attempted . 

Experimental Program 

The nature of the present investigation immediately introduces the 
concept of conditional sampling, as a way of evoking the particular portions 
of the flow of most interest. A review of the techniques available, as well 
as the specific method to be used in the present experiments is outlined 
below. However, it is important that the behavior of the jet flow in the 
mixing and self-preserving regions based on available experimental data, be 
investigated; in particular, the fluid dynamic concepts associated with 
the large scale structures. This not only helps in the understanding of 



the various nechanlsns Involved but also places in perspective the problems 
that may be encountered. 

Speculations about the existence of large scale structures in free shear 
layers developed from the two-dimensional experiments of Brown and Roshko [A]. 
While investigating density effects, they observed an organized coherent 
structure, shown on shadowgraphs, persisting downstream of the splitter 
phase at a Reynolds number of order 1 x 10^, Winant and Browand [5] 
earlier showed that the vortex pairing process in a two-dimensional shear 
layer was a result of shear layer instability, but the Reynolds number in 
their experiments was quite low. Later, Dimotakis and Brown [d] investigated 

g 

the existence of large structures at a Reynolds number of 3 x 10 . They 
verified the existence of the previously observed structures and showed that 
there exists a periodicity to the phenomenon. However, they also showed 
the development of three-dimensional eddy structures and concluded that 
initial conditions as well as the finite extent of the apparatus play an 
Important role. Chandrsuda et al. [7) also pointed out the three dimensional 
development of these structures and their sensitivity to upstream free stream 
turbulence. 

The interesting outcome of the above work is the applicability of these 
ideas to noise generation in the mixing region of a jet. The question that 
has been posed is, whether the large scale structures are direct or indirect 
noise generators? In the latter case, how does the excitation of large 
structures in a jet result in increases in broadband jet noise radiation? 

The above questions are fundamental, therefore, in guiding the experimental 
procedures. Considerable insight, closely associated with our program, is 
provided from the experiments of Yule [8]. He studied the structure of the 



mixing r<diglon at a Reynolds number of 2 x 10'^ and concluded that the observed 
vortex rings are only a transitional phenomenon and also that there ave 
substantial differences between these vortex structures and the structure 
in the fully developed region. He established that the coherent eddies are 
dominating components but are far from deterministic. Although a wide range 
of data are available on the acoustic characteristics of excited jets, very 
little has been reported on the flow characteristics. A thorough explana- 
tion of the phenomenon will only be provided* however, from measurements 
inside the flow itself. Hussain et al. [9,10,11] have reported flow 
measurements, but so far they have concentrated on low excitation levels and 
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Reynolds numbers of up to approximately 2 x 10 , while Klbens et al. [12,13, 
lA] have concentrated on the behavior of the vortex pairing process, 
which appears to only occur at very low Reynolds numbers. 

In our program, therefore, a major effort is being made to provide 
(1) detailed turbulence characteristics under excitation, (ii) a comparison 
with unexcited tests and (iii) a comparison with analytical models which 
is presently lacking. 

As outlined above, a powerful tool for detecting specific flow events 
is conditional sampling. The use of this technique presupposes, however, 
certain criteria, and the choice of these criteria is a fundamental and 
critical process. Yule [8] after obtaining velocity time histories, chose 
peaks in the velocity fluctuations as the criterion. Velocity levels were 
also used in the experiments of Bruun [15], In excited jets, sampling can 
be initiated from the forcing mechanism providing the trigger (phase- 
averaging) or r, choice of certain level of intensity at the exit of a jet 



being a fraction cf the measured non-exclted value [9], A novel and power- 
ful pattern recognition technique has been developed by Wallace et al. [lb, 

17] • It has been primarily applied to boundary layer flows and Is a digital 
method utilizing a fully computerized system. Short time temporal averages 
are selected and criteria are based on sharp gradients of velocity and 
acceleration detected from digitized output signal histories. Blackwelder [18] 
and Wygnaskl [19] also outlined concepts In pattern recognition of coherent 
eddies. 

In the present program a number of the above Ideas will be Incorporated 
based on an overall philosophy as shown on Figure 1. 

The following measurements and processing techniques will be used: 

1. Velocity time histo’Xes, longitudinal and transverse under no 
excitation using hot wire sensors. 

2. Same as (1) but under excitation. 

3. Pressure signals in the near field. 

A. Digitization of signals under (1), (2), and (3). 

5. Search for clearly defined events such as peaks or other 
repeatable values. 

6. Selection of the conditional sampling criteria and comparison with 
other established criteria detected under the same procedures above. 

7. Turbulence characteristics and phase averaging of first and second 
order moments of the velocity. 

8. Comparison of excited and unexcited cases enabling deduction of 
'Vave induced stresses." 

9. Comparison with theory. 



The experimental facilities are nearing completion* The circular pipe- 

nozzle system has been installed as well as the traversing gear. Automatic 

positioning of two coordinate axes is possible, namely transverse and 

vertical as shown in Figure 2. For the transverse motion (x) 10,000 step- 

motor steps are provided per 1/2 in. and for the vertical (y) direction, 
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10,000 per in. Accuracy can therefore be obtained within 5 x 10 inches 
for X and 10 inches for the y direction. The motors are driven by SLO-SYN 
Indexers, and the downstream positioning (z) is performed manually. 

Manufacture of the exponential horns is complete and the exact 
dimensions are shown in Figure 3. Delivery of two of the four loudspeaker 
drivers is still pending and experiments will be conducted using only the 
system of two speakers. The processing instrumentation is shown in Figures A, 
5, and 6. During preliminary tests to establish the proper operation of the 
installed nozzle system, an exit velocity of 210 ft/sec (64 m/s) maximum 

5 

was measured, yielding a flow Reynolds number of approximately 4 x 10 . 

This value is lower than anticipated and reasons for this can be 
attributed to heavy losses in fan efficiency at high speeds. Considerable 
Interference has also been encountered close to the operating facility due 
to other experimental rigs in the vicinity of our own. Heavy vibration 
loss dominated the interference resulting in contamination of the hot wire 
signals. Recirculating flow has also caused problems. The facility will 
be operated during evening hours to eliminate the external influences 
mentioned above. A program has been Initiated for the construction and 
establishment of the jet facility proposed and shown in Figure 7 , 

It is expected that flow measurements of turbulence properties of the 
jet in unexcited and plane-mode excitation will be completed by September 
1980 and that phase averaged measurements for the same excitation conditions 
will be available by the end of the grant year. 
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FIGURE 4 : Instrumentation for Fluid Dynamic Measurements 
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Abatract 

A aodel la propoaed Cor the change in turbulent 
atrueture of a round Jet In the preaence of an 
aeouatic excitation. The excitation la aaaumed to 
trigger inatabllicy wavee of a known Initial 
anplltude at the Jet exit. As theae vavea propagate 
downatreaa they extract energy from the mean flow 
and transfer it to the randoa turbulence. This 
results in an Increase In the levels of the turbu- 
lence and a resulting Increase In the radiated 
broadband noise. No calculations are presented 
for the noise radiation • however > an exanlnatlon 
is asde of the effect of excitation level and 
frequency on the Jet Clow. The nuawrical procedure 
allows for radial as well as axial variations In 
the averaged properties of Jet to be calculated. 

The results indicate that the presence of a finite 
aaplltude instability wave increases the spreading 
of the Jet. It does not vary the characteristic 
radial shapes of both the axial mean velocity and 
the turbulent kinetic energy. An energy budget 
for the tandom turbulence shows that it is fed 
energy from Che excited wave predominantly on each 
side of Che Jet lip line. This results in a broader 
radial shape for the turbulent kinetic energy. 


Introduction 

It Is now generally acknowledged that turbulent 
shear flows, particularly free shear Clows, possess 
a large scale structure. However, the role that 
this structure plays in the generation of sound by 
turbulence and even the nature of these large scale 
nxitlons such as their orderliness, particularly at 
high Reynolds numbers, remains an open question. 

It appears that at high subsonic and supersonic 
velocities In Jets and mixing layers the large scale 
motion may be associated with wavelike Instabilities 
of the primary Clow. Since Che phase velocity of 
these travelling waves Is of Che order of Che speed 
of sound in the surrounding medium the Instability ’ 
waves radiate sound, of a highly directional nature, 
very efficiently. Calculations by Tam and Morris^ 
and Norris and Tam^ show that successful predictions 
of radiated noise from supersonic jets and mixing 
layers, at angles In the far field that Include the 
peak polar noise angle, may be made assuming that 
the radiation comes from Instability waves. When 
the characteristic mean velocity Is subsonic the 
role of the large scale structure In noise radia- 
tion is less clear. At these velocities the phase 
velocity of instability waves. If such waves are 
a reasonable model of Che large structure. Is much 
less chan Che asdilent speed of sound, so that even 
allowing for spatial variations In the wave ampli- 
tude and phase velocity (see refs. 1 and 2) these 
waves are very inefficient radiators of sound. It 
appears more likely that It Is the degree of 
disorder llness of the large structures and the fine 
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scale turbulence which controls subsonic noise 
radiation. 

However at subsonic, and most likely supersonic, 
flow velocities the large acale structure can effect 
noise radiation, In an indirect way. If large scale 
motions are excited in the Jet by, for example, 
acoustic excitation at the Jet exit then the sound 
radiation of Che Jet increases. Surprisingly, even 
if the excitation is at a slngl«< frequency the 
Increase In noise radiation tu.y occur at all 
frequencies. This phenomenon has been observed by 
Bechert and Ffitenmsler^, Mooted and Schmidt^. 

Moore, Schmidt, and Vlasov and Glnevskly” also observed 
that the levels of turbulence In these excited 
Jets also changed. The variation in the turbulence 
levels and the radiated noise was found to be a 
function of excitation level, frequency and 
aclmuthal mode number. At high excitation 
frequencies there was some reduction in the broad- 
band noise radiation. 

In this paper a model is presented for Che 
mechenisms by which the structure of the turbulence 
Is roodifiel by a pure tone excitation. The basic 
formulation Is similar to those used previously by 
this auth<ir and others, most closely to Chan^ to 
study the characteristics of wave-llke disturbances 
In turbuliVnt Jets. Howe','er this paper Is concerned 
with the variation of the random component of the 
flow. The presence of a periodic component of 
finite amplitude Is able to modify both the levels 
of the random turbulence and the mean flow charac- 
teristics. Previous analyses, because of the 
obvious simplification introduced, have made use of 
integral forms of the momentum and energy equations 
for the mean, periodic and random turbulent motion. 
Though this has permitted the axial variations In 
Integral flow properties to be calculated radial 
variations are suppressed. For example the mean 
velocity profile Is represented by a shape function 
which remains independent of the excitation level, 
etc. The present formulation permits the radial 
variations, if any, to be computed. 

Though the analysis of this paper is more general 
than previous work there remain closure psroblems 
and, in the absence of the necessary guidance yet 
to be provided by experiments, some simple schemes 
have been adopted. Thus this paper also provides 
a framework within which Improved closure models 
may be tested. 

The model to be explored in the present paper 
argues chat the acoustic excitation stimulates 
instability waves on the jet column. The amplitude 
of these waves will be determined, initially, by 
the level of the excitation. Tam°*^ has developed 
an analysis to determine the amplitude of the 
excited waves in the linear region, where the 
amplitude of the excited wave is linearly propor- 
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clonal Co Che an^llCude of the cxclCaClon. Tito 
coupling Chat occura between the acouatic excita- 
tion occura over Che entire region adjacent Co the 
nixing layer (all along the potential core in the 
Jet flow for internal excitation). However because 
of the Bismatch in phase velocity between the 
Inatability and sound waves the coupling ,ia greatest 
close to the jet exit. Xn the present analysis the 
full coupling procedure as proposed by Tam will not 
be used, rather the level of the instability wave 
at Che jet exit will be assumed to be known. There 
is no reason why Tam's analysis could not be incor- 
porated: however, because of the iterative nature 
of the present aolution and the axial variation in 
the shape of Che mean velocity profile it would be 
computationally time-consuming. 


The presence of an instability wave of finite 
BUgniCude performs two roles. Firstly it extracts 
energy from the mean flow in the initial region of 
Che Jet though some is returned as the wave decays* 
Secondly it interacts with the random turbulence. 

This interaction involves both the generation of 
additional turbulent energy and its transport. The 
increase in the level of the random turbulent kinetic 
energy causes a more rapid spreading of the Jet flow 
through an increase in turbulent stresses and, it 
can be argued, an increase in the broadband noise 
radiation. No details are provided in the present 
paper of the anticipated Increase in the broadband 
noise radiation. The paper is more concerned with 
changes in Che turbulent structure. If the noise 
producing volume of the Jet is relatively unchanged 
by the presence of the excited waves then the 
increase in turbulence levels is likely, using 
scaling arguments, to yield increases in the 
radiated noise. Once the details of the modifica- 
tion to the flow are known the noise changes may 
be calculated. However, If no details of the flow 
changes are given then no estimate of the noise 
variation is possible. 
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where an overbar denotes a tinc-lhdc'pendent component, 
a prime denotes a random turbulent time-dependent 
component and a tilde denotes a periodic component 
which is related in phaae to the acoustic excitation 
at the nosrle exit. The time average of a variable 
ia defined as 

T 


him 4 / 

T<*00 O 


f(t)dt 


and the phase average is defined by 


1 " 

<f> - Lim ^ X f(t + nt) 


N-w ^ n-1 


(«) 


(5) 


where T is the period of the excitation. It is 
readily shown thst all components of u* satisfy the 
same continuity equation, 

“»i “ “’i “ “'i " ° 

Taking the time average of equation 1 and using 
equation 6 gives. 


J P 9^3 3 3 

(7) 

where the Reynolds number is assumed sufficiently 
'high for the viseeuB Bttesses to be neglected. The 
* random disturbance momentum equation is obtained by 
taking the phase average of equation 1 from the full 
equation, giving, 


3u 


8t 


- - u'^ujj -(u'^u'^j - <u'^u’i^j>) (8) 


Mathematical Derivation 

Basic Equations 


The equation for the turbulent kinetic energy is 
obtained by multiplying equation 8 by ginu'"’ 
averaging. After some manipulation the resulting 
equation is 


In the formulation of the present analysis several 
simplifying assumptions have been made; these are 
described at the appropriate point in the text. It 
is useful to start using generalized coordinates and 
then carry out the simplifications provided by the 
geometry and the physical nature of the problem at 
a later stage. The incompressible momentum and 
continuity equations in tensor form are 
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( 1 ) 


and . 

- 0 ( 2 ) 

'here u^ is the r.ontravariant velocity tensor, p is 
tne pressure, g^J is the metric tensor and the 
usual notation has been used for covariant differen- 
tiation, The instantaneous velocity and pressure 
are separated into three components 

- u^ + + u*^ , (3a) 


and _ 

p » p + p + p’ (3b) 


u^k.j 


iIQ • j 

- u’ u’ 


m,J 


- - (u'^p') 

P ‘ .11 


- (u'\’) 
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where k' * 1/2 u'^ul_is the fluctuating turbulent 
kinetic energy, k “^k* and e is the viscous dissipa- 
tion term, f? is the "wave-induced stress" given by, 

pf^ ■ g^j^p{<u'*^u'^> - u'*^u'^} (10) 


It represents the difference between the Reynolds 
stresses, pu'^u'^ in the presence of the periodic 
disturbance and without' it, since performing a 
phase average provides the sum of a time independent 
component and a periodic component. The first term 
on the right hand side of equation 9 represents the 
production of turbulent kinetic energy from the mean 
flow, the second and third terms represent turbulent 
transport, the fourth term is the interaction 
between the random and periodic fields (this is not 
a transport term and will be shown to always 
represent a net gain in turbulent kinetic energy), 
the fifth term is the transport of turbulent kinetic 
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m 1 — 

where q ■ 7 la tha wave klnatic energy, q ■ q , 
and the vlacoua taraa have bean neglected. 


Equationa 6, 7, 9 and 11 provide the franework for 
tha praaant analyala. In the next aection the 
gaonetric ainpllficationa and the turbulence cloaure 
achanaa are introduced. 


where 1 ia teken an the width^of the jet niidng 
> layer between radii at which u ■ 0.9 and 0,1 
raapectively. A gradient tranaport hypotheala ia 
uacd for the turbulent diffunion, 

-v’(p’ + k^) • ^ • 06) 


and dlnenaional argunenta provide the energy dlaaipa- 
tlon rate aa 



(17) 


Turbulence Modela 

Cylindrical polar coordinatea are uaed with the 
origin at tha center of the circular jet noaale and 
the a axle aligned with the jet centerline. The 
phyaical coaponenta of the contravariant tennor 
are v, w/r, and u with reapect to the r, 4* *nd z 
coordinatea. After aaveral alapllficationa based on 
the ansuaptiona that the tine-averaged flow is 
axisynawtrlc and axial variations of tine-averaged 
values are negligible with respect to radial varia- 
tions the following systca of equations is )btalned. 


- 3u . - 3u 1 3u*v ■ r 1 3aVr 
“ 3z 3r ” r 3r " r 3r ' 


( 12 ) 




v^(p' + k’)r 


u'v' ^ - c 


where 




(13b) 


The integral form of equation 11 which will be used 
in the subsequent analysis is given by 

00 ® ® 

/ uqrdr “ “ / “v rdr + / rdr . (14) 

o o o 

In Che above equations the velocities have been non- 
dlnenslonallzed with respect to the average jet exit 
velocity u., lengths with respect to the jet radius 

2 

rj, and the pressure with respect to PjUj where Pj 

is the jet exit density (equal to the density 
everywhere in the present formulation) . 

Examination of the equations above Indicate Chat 
several turbulence models are required. In 
equation 12 both the random Reynolds stresses and 
the averaged wave stresses must be modelled. In 
equation 13a the turbulent diffusion, viscous dissi- 
pation, wave Induced stresses and "wave transport" 
must be modelled. In the present analysis It is the 
modelling of the latter two effects which are of most 
Interest thus a conventional "one-equation model" 


The coefficients C^, C2 and Oj^ take the values 0.05, 
1.50 and 0,7 respectively.. These models would 
enable Che jet development to be calculated in Che 
absence of an excitation. Row consider the "wave- 
induced stresses." By analogy with the Reynolds 
stress model they could be written 

rJ ■ - 2 S^sJ , (18) 

where s|^ is Che periodic rate of strain tensor, 

“1 ■ 2^“’i 

The choice of ^ will be discussed below. Contrac- 
tion of equation 18 leads Co 

r^ - 2(<k'> - k) - -2 c^sj . (20) 

However from the continuity equation, s^ ■ 0, so that 
’’ the random kinetic energy In the presence of an 
.rjicitatlon is equal to its value without excitation, 
which Is in disagreement with observations, Thus 
equation 18 should be written, 

r^ - -2 e^sj +|(<k'> - k)6^ . (21) 

Now consider the choice of Ct-. The value of C Is 
given by equation 15b as k^'2 j;,, 
consequence of the assumption that there is only one 
characteristic time scale for the random and mean 
flows (see Tennekes and Luralcyl®) . Letting this 
time scale be t the eddy viscosity may' be written 

\ - Cj kt (22) 

If Che time scale for the periodic field Is denoted 
by t the periodic eddy viscosity could be written as 

\ ■ ^4 » C23a) 

or 

CjE^/Et - t/t , (23b) 

where C3, Cl and C5 are unknown coefficients. Now 
ujt/rj ■ 2n7(a ■ 1/f, where is end f are the non- 
dimensional radian frequency and frequency of the 
periodic excitation, respectively, and t £./k^/2. 

So that 
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i 


. 1/2 


B. 

2t 


( 24 ) 


where! the Strouhel nui^er, St ■ f*dj/Uf, where f* 
ie the diacneional excitation freuqency and dj 1 h 
the jet dlaaeter. ThuBp 

\ " 2^ • 


Since the value of k varlea little with axial 
dletance the periodic eddy vlecoalty la nearly 
conetant for any Riven frequency, and decreaaea with 
Increaalng frequency. 


In the preaent paper the "wave-induced atreeecB" 
are aodelled ualng equation 21. However the 'Vave 
tranaport" term In equation 13a la neglected, Thin 
tern would have been aero If the aodel equation IB 
had been weed . Thl a neana that a further model for 
the quantity v<k’> la not required. It should be 
noted that Chan/ uaed equation IB with the rusultinp, 
neglect of wave tranaport. However in hia analysie, 
olnce he did not dlatlngulah between aource and 
tranaport teraui, the wave Interaction terms In 
hia analyala contain many termB which would be 
eliminated using the continuity equation. In the 
present paper calculations have been performed 
asaumlng that \ 


From equations 13b, 18, and 19 an expression for 
^ may be obtained. In the present p/iper only axl- 
syametrie periodic disturbances have been considered. 
In this case alapliflea to. 


^ i 4 2 


8vi 2 . 3y2 j 


( 26 ) 


Note chat 4> la always less ctun zero so that it 
represents a gain of turbulent kinetic energy and 
a loss of wave kinetic energy. 

In light of the dynamic Instability of the jet 
mean velocity profile It is reasonable to model the 
periodic flow field as large scale instability waves 
stimulated by the excitation at the Jet exit. These 
Will be modelled as linear waves whose wavelength 
and growth rates are the eigonvalue of a local 
invlscld stability calculation. So that a periodic 
disturbance la written, for example as, 

u - A(z)Re{Q(r) cxpfKaz - ict)]} , (27) 


where a Is the axial wavenumber and A(z) represents 
an axial amplitude variation. This can lead to some 
difficulties, particularly when the wave is decaying. 
In such cases the invlscld solution is not valid 
over the entire real radial axis. (Physically 
correct damped invlscld solutions do exist. However 
they require a contour deformation Into the complex 
r-plane; see Tam and Morris^.) In order to over- 
come this difficulty the instability wave character- 
istics will be obtained from a local viscous 
stability’ calculation. The Reynolds number used In 
the calculations is suf flcleutly high that the solu- 
tions apnroach the invlscld limit. Following Tam 
and Chen^i the value was chosen to be 500. Since 
It is the intent of this and subsequent analyses to 
consider both incompressible and compressible flows 
(to examine amplification effects In high Mach 
number Jets) a technique has been 4eveloped to 
consider the stability of viscous, compressible 


•iF runli quality 

axiaymmctrlc jets, Tliough only results for axi- 
symmetric disturbances in the incompressible limit 
are given in the main body of the paper, the 
complete analysis is given in Appendix A. 

Tlie stability snslyals gives all the unknown 
quantities in equation 27 with the exception of the 
amplitude, A(c). This is obtained from equation 14 
which now representa an ordinary differential 
aquation for A(x) . 


1*1 


h>'-h-h 


where 


1 t - 


/ u(lu| + |0| )rdr , 


I2 • ^ / Re{CiO*}-~ rdr , and 


(28) 

(29a) 

(29b) 

(29c) 


where an asterisk denotes the complex conjugate. In 
order to provide a consistent definition of A(z) 
the eigcnsolutlons u, v and p arc nonnallzed such 
that 

- 1 at r - 0 . (30) 


Numerical Solution 

The solution of equations 12, 13a and 28 has 
been obtained numerically using a technique which 
employs a standard three-point, variable step-size 
finite difference formulation in the radial direc- 
tion with the axial derivative discretized in a 
manner which can either be fully Implicit or 
explicit. It is hoped that sufficient details of 
the numerical scheme are provided below to enable 
the reader to use the technique. (It has been this 
author's experience that the distance between a 
numerical scheme and a working program is a very ^ 
long one.) 


Introducing the turbulence models given in Che 
last section the axial momentum and turbulent 
kinetic energy equations may be written, 


- 8u 
“ 3z 


+ v 


3u ^ ^ 3uvr 

3r ^ r 3i’^ ^ ^ 3r'’ ” r 3r 


(31) 



(32) 


12 

Drawing on the work of Spalding and Patankar and 
Chan/ new coordinates (C,w) arc Introduced such 
that 

C - and <i) - ~j7^u(z,s)BdsJ^'^^. (33) 


2 

he is a measure of the total axial volume flux In 
the Jet. Due to entrainment this Increases with 
axial distance and must be computed at each axial 


4 



, location t The celculotlon of d(n|)/d> !• deecrlbed 
below, in tcrM of thee# new coordinetee equetlona 
31 end 32 bccone, 


IS 

H 


IS 

3C 





< 34 ) 


2m ‘‘<1^ » 

^2 dc 8u) 


+ i, it -H. isi 

(i» 3u‘ u) Ok. 

oTT^ ” "T* 


+4 


u .3ut 7 - 


(35) 


PAGE IS 
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reMining aetrlcee oc^ given In Appendix B. The 
boundary conditions ace that the axial velocity and 
turbulent kinetic energy are aywMtric about r ■ 0 
and they both vanish at large r. This leads to valuee 
Si* ^1* di i ■ 0 and N which are also 

given in Appendix B. 


Equation 40 and the additional equations provided 
by the boundary conditions uy be solved by a 
standard algoritha for inversion of a block tri- 
diagonal utrix. However the axial variation of 
(hg) and the initial conditiona Bust first be 
specified. 

2 

Calculation of d(ri.»)/dg 


The flow field is divided into lines of constant 
hi denoted by the subocript i, 0 ^ i d N, and lines 
of constant ; denoted by the euperscrlpt j. Then 
at any point (i,J) the equations sksy be written. 


i3u\^ 

'3?^i 


a2„ 

{a-^ + b 


and. 


3h)‘ 


Is 

3h) 




^35^1 


( a 3^k 
\ 30)^ 


+ d f + ek + f } 


j 


(36) 

(37) 


1 1 

The definitions of the functions aj through f^ ate 

given in Appendix B. The derivatives with respect 
to u are written in a finite difference form such 
that 


—t; 2 (®i+l “ ■*■‘^“1-1^ 

. 3o> 2 1 (l+a)a2^"^(h,)^ ^ 

* /QB 


and 


(|^) ^ i (*i+l * ° ®i-i) 

i (1 +o)aS^ ^ 


(38a) 


(38b) 


where s is cither u or k, h. is the stepsize given 
by (h)^ - ‘^l-l^ ° the^step-size ratio h^/h^_j^, 

which is constant. Setting. 

the axial derivative is discretized in the usual 
finite-difference fashion. 

" ®i “ 

where Ai; is the axial atep-size. If X ■ 0 this 
is an explicit form. X ■ 1 is an Implicit form and 
X - 1/2 is the Crank-Nlcolson form. The last was 
used in the present calculations and was found to 
give rapid convergence. Substitution of these 
forms into equations 36 and 37 enables them to be 
written as either a pair of tridiagonal matrices or 
as a single 2x2 block trldlagonal matrix. The 
latter form, used in the present work is, 

-i -i-1 -i -i^^ -i -i+1 “ -1 ’ 1 -= 1, N - 1 

(40) 

! . - T 

where the solution vector ■ (u^, k^) , and the 


At the ouf»r edge of thn Jet 3u/3C 0 and w 1. 

Thua from equation 34. 



(41) 


Following Spalding and Fatankar^the second term is 
Ignored and the last term will be allowed for by 
writing, 

d(nf) a 

--2f--Lim {(1 -ke) (Ej.r)} (42) 


The right hand side is approximated by its value 
at ti) ■ Then if. 


G(ai,z) - -(1 + G)ej.{l + (43) 

n.. 


Then, 


‘Xnf) 

dz 


a - wcj., + »5J« 


(44) 


and. 

(Dg)^’*’^ - (n^)^ -k 4A?(1 - X)gJ_^ + 4a?xgJ+J , 


(45) 

where the tilde denotes the estimate of G from a 
previous Iteration. 


The Initial Conditions 


The mean velocity profile at the jet exit was 
taken to be the developed mixing layer form given 
by Maestrello and McDaidT-^, 


u(r,0) - 


exp [-15. 3825 (r - 0,796)^) 

1.0 


r > 0.796 

r < 0.796 
(46) 


Choice of this developed profile for the shape at 
z - 0 does ignore the short region of adjustment 
at the jet exit from the boundary-layer-like profile 
of the Jet nozzle. The initial value of H is then 
0.304. The turbulent kinetic energy profiles were 
estimated using a constant eddy viscosity assumption 
which leads to 
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• k - 0.5!277(t - 0.796)u , r > 0.796 (47) 

In thii potantlal corn rnglon th« tucbultnt klMtic 
•nnrgy wmi n«t to a conntant valua of 6 x 10'^ which 
corraaponda to a typical axit turhulanca intanaity 
of 2 parcant. 

Thaae Initial conditiona though raadily apecifiad 
) ara aoMwhat idaaliaad particularly in the potential 
cora and reault in a non-raaliatic behavior of the 
jet in Che potential core region. Fron equation 35 
it can be aaen that if both k and Q arc conatant 
the turbulcnea can only decay. Thia waa obaerved 
in the calculation and reaulted in an cxtenalon of 
the potential core length. However thia paper ia 
concerned with nodificatlon to a baaic undiatuebed 
flow by an acouatic excitation rather than abaoluce 
behavior thua no effort haa been made to correct 
thia behavior. Such a correction will be made in 
later calculatlona. 


With u apecifled the variation of w with r may 
be calculated from equation 33. The valuea of u 
at the grid valuea of u were c)ien determined uaing 
a cubic apline fit for Q(u). 

Calculation of t fron u and u(ie ) 

At each axial location the valw of radial 
poaitlon nuat be calculated fron the local valuea 
of 0) and G. Letting f (u) ■ u/u and using 
equation 33 

- (r^)^_j^ + 4n^hi / f(t)dt , (48a) 


where 


t - (w - • 


(48b) 


A cubic polynonial spline fit may be obtained for 
f(t) on the interval <*>£) ®nd finally, 

2 2 2 ^1 
(r^)^ - (r^)^^l-t-4n;h, f 


(49) 


where 




(50) 




The Computer Proeram 

The nunerlcal method was programmed in FORTRAN 
and run on an IBM 3033 processor. A simplified 
flow chart is shown in Figure 1. After initializing 
the values of u and k at the grid points the elgen- 
sdlutions of the Orr-Sommerfcld equation are 
obtained for that velocity profile. This gives 
the distributions Of the periodic flow field. The 
integrals in equation 29 may be calculated as well 
as the radial distributions of the source term 
and the periodic stress gradient (3(uvr)/3r]/r. 

The wave amplitude at the next location is then 
computed from equation 28. Estimates of the 
coefficients of the matrices in equation 40 may 
now be made and the values of u^"^^ and 

calculated using the matrix solver. The wave 
sha^r. and Integrals at the downstream location 
ara then calculated from the Orr-Sonnerfeld equa- 


ap PUMit (QUALITY 

tlon using the ostlmated downateaxm velocity 
orofilo. Reviaad •atimates of tha wave amplitude 
my also be mado and new estimataa of the values 
of the matrix elfementa in equation 40 are made. 

The procedure la repeated until the mean velocity 
profile at the downstream location la unchanged by 
further iterations. This requirement was met in 
the present calculations when the aun of the root 
mean square differences between the new and old 
estimates of at all the grid points was less 
than 0.011. Thin was found to ensure three decimal 
places of accuracy at each grid point. 

The numerical method ensures that alterations 
in the mean flow properties due to the presence of 
a finite amplitude periodic structure can be 
observed in both the axial and radial directions, 
and that the "shape” of the mean velocity and 
turbulent kinetic energy profiles may alter. The 
grid constants and coefflciente used in the present 
calculations arc shown in Table 1 where 

m a 4 ?^ . 
z 


Table 1 Constants used in computer program. 
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Fig. 1 Flow chart for computer program 
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RaBiiltit and Dlwcuawton 

Sone pr«Unin«ry calculationii h«v« b««n parCoriMd 
uslnfi tha Bodal daHCclbcd in th« previous •«ction. 
Th»»e calcuUtiona «r« vitw«d «» «n Initinl teat o£ 
the nethod to provide guidance ae to inproved BodolSt 

The axial variation of the Jet thickneee Cot an 
unexcited Jet ia ahown in Fig. 2, T)i« variation 
approachea a linear increaae with axial diatance 
with a virtual origin of x * -2.0. Title aereee with 
the neaaucenenta of Maeatrello and McDaid^J whoee 
data were uaed ae the initial velocity* proCile, 

A collapse of the conputed Man velocity and turbu- 
lent Kinetic energy proCilee la obtained by using 
k Btretched radial coordinate, (r -1)/(X'1‘2) . 

Tlieee are shown in Figs. 3 and 4 respectively. Thus 
the unexcited Jet behaves in agreenenc, with experi- 
■ent, though the spreading rate shown in Fig. 2 
is somewhat low due to the idealieed nature o£ the 
potential core mean velocity and turbulent kinetic 
energy profiles. 

if 

Moat of Che calculations hav>.’ been perforsied for 
an excitation Strouhal number of 0,25. However in 
Fig. 5 the effect of Strouhal number of the Jet 
spreading is shown. The spreading is greatest 
for St - 0.5 and least for St • However the 

excitation level of |a|q " 1 x 10**^ was the minimum 
level at which modifications to the Jet structure 
were observed. In all cases the Jet i^read rate 
increased. Fig, 6 shows the effect of increasing 
the excitation level by n factor of 10. The Jet 
spreads initially more rapidly and then grows 
linearly for an excitation Strouhal number of 0-25. 






Fig. 5 ■effect of excifgeion fieequoncy en joi 

thickneiB. |a1q "lx 10"^. St "0.1,— — ; 
St ■0.25,-——; St “0.50, — — — . No 
excitation * 



Fig. 6 Effect of excitation amplitude on jet 

thickness. St ■0.25. |Ap “lx 10"^, — 
[A|q ■lx 10"3, — . °No excitation, 



Fig. 7 Variation of relative amplitude of wave, 
iAr/|A[", with axial distance. St ■0.25. 

?or legend see Fig, 6 . 

The new virtual origin of mixing is at x ■ "3»75. 

' The absolute levcla of excitation needed to produce 
these effects are discussed below, 

When there is no interaction between ‘the wave 
and the flow field the local rate of growth is 
that predicted by invlscid linear stability theory. 
However, when the wave is of sufficient amplitude 
to generate significant 'Wve-induced stresses'* the 
wave loses energy to the random turbulence. This 
can be seen in Fig. 7 where the ratio of the wave 
amplitude to its initial amplitude is shown. From 
the definition giving |Ap using the normalization 
of equation 30 it can be associated with the mean 
square pressure fluctuation on the jet axis. As 
the wave amplitude increases its rate of growth 
falls below that of the linear prediction. For 
an Init'i.al amplitude of 1 x 10“3 for x/rj greater 
than 4.0 the wave transfers more energy to the 
random turbulence than it gains from the mean flow 
and it begins to decay. There is also some decrease 
in the growth rate due to changeo in the shape of 
the mean velocity profile and the jet width. How- 
ever for 1 a|^ ■ 1 X 10“^ linear stability predicts 
wave growth up to x/rj ■ 5,5. 

The amount and distribution of the energy trans- 
ferred to the turbulence is seen by performing a 
radial energy budget. This is shown in Fig. 8 for 
the unexcited case, for a jet width of 0.447 which 
occurs at x/rj * 1.66, Fig. 9 shows the energy 
budget in the excited case at a similar jet width, 
0,446, which occurs closer to the jet exit, x/rj"1.07, 
in this case. The additional term in the turbulent 
kinetic energy equation which describes the Wave 
interaction effect is seen to peak on each side of 
the jet lip line at r/rj ■0.8 and 1,2, Its peak 






Iil2t 




B ' ‘ ‘ ^ - » t j 

• ,S 9,6 t.7 t.l t.t l.t l.l 1.2 1.3 


Radial dlatanca v/tj 

rtg.‘'8 Enargy budgcL i'-' ' no agcltatlon. 

width ■'Oi447. it'fi.j ■l.fifi. Eroductl^n, — 
Dlfluaion^ — — S Dlaaipation, 

Wave Interaction^ — — , 



Radial diatance r/tj 

Fig. 9 Energy budget with excitation. St*0.25| 
|a(” "lx 10“3, width ■O.A46, x/tj-1.07. 
For°legend nee Fig. 8. 
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value ia an order of »iagnB^<Plt^l^tfty4lleIPy 
direct production Iron the Mean flDU» however aince 
the other effecta are aliaoit in equillbrlua* moat 
of the production balanced by viacoua dieelpattoni 
thia haa a aignificant effect. 

the mean velocity profilaa in the excited caae 
are ahown in Fig. 10. T)ie ahape of the wave profilaa 
In aueh chat in the potential core the wave tranafera 
'Monentun to the aean flow. Thia reaulte in a alight 
increaae in the jet centerline velocity. The data 
all collapae uaing the aaae atratehed coordinacci 
(r -ll/Cx-fT), aa weed for the unaxoited jet, with 
the exception of the profile at x/rj"2.0* How> 
ever, thia ia to be anticipated aince in thia 
region the jet growth rate ia nore rapid, An 
initial look at the velocity ahape, in coapariaon 
to Fig. 3 ahowa a broadening of the profile at the 
outer edge of the jet in agreeoent with the Mcaaure- 
■enta of Schnidt.^ However if a new virtual origin 
ia choaen at x/tj ■ 3,75 the calculationa iron both 
excited and unexcited jata are virtually identical. 
Hie aame appliea to the profilaa of the turbulent 
kinetic energy which are ahown In Fig. II aa a 
function of (r -l)/(x +3.75). The ahape of the 
profile ie unchanged Iron that of the unexcited Jet 
except at the fartheat downatrean atetione, 
x/rj“7.96 and 10. 0. In theae locationa there ia 
a rapid increaae in the turbulent kinetic energy 
ineide the lip line. At x/rj"l0.0 there ia a 
Budden increaae in the peak level to 0.033 pe«iking 
at a radial location of r/cj " 0.66. An initial 
look at the energy budget for thia atatlon indicatea 
a iudden ipae Of energy by the wave which ia elao 
evidenced by its rapid decay shown in Fig. 7. 



(r-l)/(x+2) 

, Fig. 10 Collapse of axial mean velocity with 

-l)/(jt +2) for excited jet. St-0.25, 
|A|q ■ 1 X 10“3, For legend see Fig. 3. 
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FiK> 11 Coll«p8a of tucbulenc kinetic energy with 
fr -IV (x+ 3. 75) ^ cx-eitea jet. St ■ 0 . 25, 
lA|‘ - I X 10“3. For legend see Fig. 3. 



Fig, 12 Effect of excitation on radially integrated 
turbulent kinetic energy. For legend see 
Fig. 6. 


Hie effect of exeltation elec increaeee the 
integreteo turbulent kinetic energy given by, 

m 

/ krdr * 

0 

Ttile ie etiown in l^igt 12 where It can also be aeon 
chat at tha lower excitation level there la e 
ellght decreeeo in the integreCkU levels. 

The choice of excitetlon level in theae calcula- 
tions hsa baen daesnalnad by whatever value gave 
appreciable changes, Aa notad above the nor«atiaa- 
tion i»ed for the elgenaolutions of the Orr- 
Soaaerfeld equation Meant that iA]| la tha Msan 
aquare treasure fluctuation at the jet exit. Thus 
for |a|o of 1 X 10***) the lower excitetlon level, 
the preseure fluctuetion et the Jet exit la 2 per- 
cent of the dynamic heed at the Jet exit (token as 
p uj/2) « Tills is nuch .i^erger than the miniimim 
level obaervad by Hoore of 0.Q8 percent. In 
order to achieve conpacabla effocta with reduced 
excitation levels it would be necesaery to increeae 
the value of 

w 

Coneluaions 

The results indicate that the model can predict 
observed effects of acoustic excitotion. However 
the absolute levels required dp not mteh exporlmencs. 
This may be attributed to the choice of model for the 
wave-induced stresees particularly its multlplicativa 
constant. Tlie results point to the need for a 
knowledge of the following: 

1. further effects at frequency and excitation 
amplitude, 

2. the effect of choice of coefficient in c^, 

3. the effect of permitting ^ to be frequency 
dependent, and 

A, the effect of Inclusion of "wave-transport" 
tenoB in the turbulent kinetic energy equation. 

■nicre also cxirts a clear need for further 
oxperl"*ntal docujiontatlon of the changes in 
turbul ; structvxe. Without this guidance, 
development of ajiproprlate modela, for wave- 
induced stresses in particular, is made difficult. 
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Appendix A 

Solution of the Viscous. Compressible, Linearized 
Stability Equations for an Axlsyametric Jet 

In the analysis It is assuaied, following Dunn 
and Lln^^ that the viscous dissipation terriB in the 
energy equation nay be neglected. Assuming that 
the mean velocity and density are functions of 
radius only the non-dimensional equations of motion 
may be written, 


90 ' . 4 . 1 + hih +v» Si -f u ■ 0 , 

(A.l) 

-/9v* - 9v'^ . 3 p1 . lfr,2„, vL.J.M'i 


n 3w' w 1 3 e 1 . JLrp2 • . slI + J. 

PllF * “ * r rI^ " r2 r^ 


<A.2) 


(A. 3) 


where the velocity fluctuationa in the r, and z 
directions are v', w‘ and u', reapectlvely, and. 


R- 

V 


where pj, Uj and rj are the jet exit density, 
velocity and radius, respectively. Also 

2 ' 


(A.5) 


31*2 t 3r pi. g^2 gp 


= {- 


(A. 6) 


Combining the energy and continuity equations leads 


to, 


„2>8pV . - MuilEi + + i » 0 , 

”j'^ + “ 3z 3z * r 3r ^ r ' ’ 

(A. 7) 

where M^ ■ u^/a"? and a. Is the speed of sound at the 
J J J J 

jet exit. In equations (A.l) through (A. 6) all 
variables have been non-dlmenslonallzed with 
respect to the corresponding Jet exit property and 
the pressure is non-dimensionalized by the dynamic 
head, p.u^, at the jet exit. Equations (A. 2, 3, 4 

and 6) provide four equations for the four 
unknown quantities, u', v', w' and p’ . 

The disturbances are written In a wave-like 
form BO that a typical disturbance, q'(r,i^,z,t), 
may be written, 


q'(r,({>,z,t) ■ Re{q(r) exp[l(az -tut -fmiji) )} 


(A. 8) 


where oi Is the axial wavenumber, tu Is the radian 
frequency and m is the azimuthal mode number. 
Substituting expressions of the form (A. 8) into 
the equations of motion gives, 

lp(otu -tu)y +p' * ^ ” [a^-f -i ^ w), 

R r pit 

(A. 9) 

lp(aS -tu)w - I |0" + ^ - (a^ +1 <f), 

r R r 

(A.IO) 

lp(au-tu)G + pu'v + lap *= -[a2+ ^]u] , 

(A.ll) 
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{ 


pnd 


lMj(aS -u))^ + {0* + i + ga(j) . 0 


iiid 


(A.12) 


ThlH expresaion reducoa to the known aolutlona In 
the Invlacld, and Incoopreaelble caaua, that la, for 


Mj-0 




(A. 18) 


and for 


Priata denote differentiation with raapect to r. 

Pov Mj ■ 0 and P ■ 1 theae cquationa raduce to the 

InconpraaBlbla cquationa. For v> 0 or R ■ the 
cquationa raduce to their Inviacid fora. A aolutlon 
for the former caac for the axlayanetric jet vaa 
obtained by Morrla^^- and for the inviacid caee by 
Hichalke^". The botaidary condltiona to be aatiaflcd 
are*' 

0, 6, Q, 0 *«' 0 aa r 

u(0) ■ p(0) • 0 n ^ 0 , 

v(0) - w(0) - 0 n ^ 1 , 

v(0) + 10(0) - 0 1 , (A. 13) 

Ae in the caae of IncompreaBible flow thla alxth- 
order ayatea of equations reduces into one fourth- 
order aystea and one second-order system for a « 0. 

That is, the w equation Is not coupled to the 
equations for Che reaaining variables. Since the 
solution in this case, a ■ 0, is most readily 
obtained the analysis is described briefly below, 
and is then used to infer the general solution 
for arbitrary values of a. 

In order to obtain a. nisBerical solution to thl@ 
problea it is necessary to deCeraine the behavior 
of the solutions in the regions of uniform mean 
velocity, that is outside the jet flow and in the 
potential core region close to the Jet axis. 

In the potential core region the non-dlaensional 
mean velocity and density are unity. A series 
solution to the reduced equations is sought close 
to Che jet axis. Thus the disturbances are written, 
making use of the boundary conditions for a ■ 0, 

V - V^r + + V^r^ + 

u - + Uj^r + + U^r^ + — 

^ * *’o V **3’"^ 

After some algebra it may be shown that, 

f , l«^R+lA^M^(«-iu)l , 

p - pJl + i r^ 

A » J ej2 J 


R » p ■ Pq^oC*^^ " (a-(ij’)^Mj • r) (A. 19) 


[Rj, + iMj(a-u))l 

j la2R+iA2„2(a_^j)j2 ^ 
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lR+lMj(ot -6))]^ 


+ 


where 


+ 1(0 - (o)R . 


(A. 15) 


(A. 16) 


This is the series expansion for the modified Bessel 
function of the first kind and zero order, that is. 



o\ + iX^M^Co -oj) 

- - - I I ^ I I 

[R + lMj(a - (ii)] 


(A.17) 


Inspection of the continuity equations then suggests 
the associated forms for Q and 0 and substitution 
of. these forms into the equations of notion yields. 


a[R + lM,(o-(A)) 

- **0 (o - A)" 


IX* [R + lMj(0 - ui)J 


<A.20) 




where 


X* 


-i' lxS!j(0l - (D) 
(R + lMj(o - «il) ] 


(A.22) 


It is clear that the disturbance prcssucc has only 
one linear solution which satisfies the boundary 
conditions. This is the so-called Inviscld solution 
for incompressible viscous analyse; bince the 
momentUB equations, which arc the i/nly equations 
effected by viscosity are identldl for both the 
incoapresslble and compressible analyses the 
"viscous solutions" should be identical in both 
cases. It is then readily shown that another 
independent solution of the disturbance equations 
is, 

p ■ 0 


6 ■ U„I„(Xr) 

« ■ -"o f 


(A. 23) 


The same analysis may be applied Co any region 
where the flow is uniform. The solutions for any 
region in which p “ and u ■ u^ may be written as: 


0: 


p - (iX*r) 


(A. 24) 


a(R-iM,n) 


u •• A 


1,2 p nR 


^ (1),(2) +A^^^H^^^^»^^\lXr) 


„ (1),(2) 


^ " “^1,2 p^nR "1 

-s.^f 


(lX*r) 


(A. 25) 


(A. 26) 


B 0: p - A 2 (iX*r) (A. 27) 

G - A, „ H H (1).(2) 

(A. 28) 
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where 






*r - I 


2^ gi+i 

(1 ^ (l+a)o"hj 


(B.8) 


(A, 29) 


* ■ *1.J « <“*') +*3,4 f 


•r‘ ■ I- “inKr — 2 *r‘ + - I 


where 


and 


(A.30) 

n ■ u > au , (A, 31) 


, 1+1 


(l+o)a^^“^(hp^ 


(B.9) 

ir'+ h-b^XAC 


(l+o)oSij 




oS - ix^i^n 

X»2 . ^ 

(R - iM?n) 

J 


(A. 32) 


(B.IO) 


(B.ll) 


Theae eolutlons have been %rrltten in a fom appro- 
priate to the conditions far from the Jet axis. The 
'’'lolce of the Hankel ftnctlone of the first or 
second kind depends on the argument and the satis- 
faction of the boundary condition at Infinity, A. 13. 
Similar forms, In terms of modified Bessel functions 
may be written for the solutions in the potential 
core or close to the J.et axle. 

Appendix B 

In equations 36 and 37 the definition of the 
functions a through f are, 

.2' 


+ (i-x)cJ}a 5 - nj 

2 

’t' ' a *a)a\ “^1“' 

(B.12) 

sJ+1 


,. 1+1 


“21::^ — 2 +sJ'^^}XAc - 1 

^ /O 1 0 


(B.13) 


~j+l 

r.1+1 ■ \ 2 

^ <I+o)a^^“^(hj^)^ \ (1+0)0^ 


^ u 

“ " t “A 

(u 


h • 1 f£_ F JL) 4. i<il 

( 1 ) 8u» '• 111 t ^ 2 


20. < 

dz ' 


2 3 


oin 


2 aul 


{uvr} • 


(B.l) 


(B.2) 


(B.3) 


k- * -TTbr 

(B.IA) 


.1 2 


-(XfJ'*’^ + (1 - X)fJ}AC - kj , 


(B.15) 


where a tilde denotes the guessed value from a 
previous Iteration. 

The boundary conditions on u and k may be wrlrfen, 


dn: 


2 E 

A i -l-f£ £. -S.) 4 . ^ 

“ 0 ) 8oiU ^2 dz ’ 


ACjk 


1/2 


u £> 


2 - 


and 




f . I__iL e 

^ 2 A ^t%/ u • 

01 Hg 

The matrices In equation AO may be written. 




u 


k - 

k. 



(B.A) 


o 

1 

o 

1 

» 

(B.16) 




- 0 

k. (K 

N 

0 



(B.5) 

Y - 
-0 

Y » 

In 

1 • lo 

- -I 

and 

!?n “ 2 • 

(B.17) 


where I is 

the Identity matrix 

and 0 

Is the 

zero 

(B.6) 

matrix7 

• 







and 

d “ 

-o 

o 

o 

v.^ 

. ^ 

- (0, 

0)’^ . 

(B.18) 


! 

X ■ i 

O ' 
< 

Y >• 

Bl 0 

-1 1 

10 A^l 

’ -1 

1 

0 B| 

5i " 

0 

and d^ - i 

(0^ , Dj) 
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Abstract 

Calculations of the viscous, spatial stability characteristics of 
a compressible axlsymmetrlc Jet have been made. The asymptotic solutions to the 
stability equations in regions of constant mean 'flow properties are developed 
and used as the boundary conditions for a numerical Integration of the 
compressible stability equations In cylindrical coordinates. Calculations 
are presented for a range of Mach and Reynolds numbers for both the axi- 
symmetrlc and helical modes. Increasing the Mach number is found to 
stabilize the flow as is decreasing the Reynolds number. The helical mode 
is found to be more unstable than the axlsymmetrlc mode for higher Mach 
numbers. The viscous calculations are' found to closely approach the 
Invlsdid limit for local Reynolds numbers greater than 500. 



Introduction 


The stability characteristics of free shear flows are of considerable 
Interest In the field of aerodynamic noise. Instability waves appear to 

perform an Important role In noise radiation at both subsonic and supersonic 

1 2 

mean flow velocities (toore and McLaughlin, et al. ) , Analytical studies 

3 

of noise radiation by Instability waves, such is that by Tam and Morris 

have assumed that because of the dynamic Instability of the mean flow 

profile in jets and shear layers an inviscid analysis may be used. Though 

this assumption appears justified by comparison' with experiment and has the 

attraction of simplicity it has a disadvantage. For damped waves, though 

valid inviscid solutions do exist, the solution is not valid over the 

entif^e transverse coordinate. If a more detailed analysis of the interaction 

between the instability wave and the other components of the fluid motion 

is desired knowledge of the wave properties at all real locations is 

required. Such interactions are thought to be responsible for the broadband 

A 1 

jet noise amplification measured by Bechert and Pfizenmaier , Moore , and 
others. Ihe analysis of the present paper overcomes this difficulty by 
using a viscous analysis. The question of the appropriate choice of 
Reynolds number in a given flow problem, whether the basic flow is laminar 
or turbulent, is not addressed. Solutions and calculations are presented 
for the viscous compressible stability equations in cylindrical coordinates. 

The inviscid stability of compressible jet flows has been studied by 
many authors including Lessen, Fox and Zien^ and Michalke^. The former 
considered a cylindrical vortex sheet and the latter used realistic mean 
flow profiles which simulated the irltial and final stages of the annular 



2 


jet nixing region. The viscouis stability of inconpresslble axisynmetric 

jet flows with realistic nean velocity profiles has been calculated by 

7 8 9 10 

Burridge', Hollendorf and Gobhart , Lessen and Singh » and Morris . All 

2 -2 

authors considered a jet velocity profile given by u ■ (1 + r ) which is 
characteristic of the flow downstream of the annular nixing region. 

Morris^^ also considered the same mean flow profiles used by Mlchalke. 

The present paper describes the viscous stability of a compressible 
axisynmetric jet to disturbances of arbitrary azimuthal mode number. The 
asymptotic solutions to the compressible stability equations are developed. 
Calculations are presented for a range of Mach numbers and Reynolds numbers 
for a mean velocity profile considered by Mlchalke^ and Morris^® which is 
characteristic of the end of the annular mixing region of the jet. 

\ 

The Stability Equations and Boundary Conditions 

The full compressible stability equations in Cartesian coordinates 

have been developed by Lees and Lin^^. A simpler set of equations which 

neglect viscous effects in the energy equation was introduced by Dunn and 
12 

Lin . A similar simplification is used in the present paper. The dissipa- 
tion terms are neglected in the energy equation and the viscous effects in 
the momentum equations are characterized by their incompressible form. 
Variations in the molecular properties of the fluid are also neglected. 
Though these assumptions restrict the analysis to moderate Mach nimbers and 
flows with little variation in the mean static temperature it should be 
remembered that the purpose of the analysis is to provide a viscous correc- 
tion to the inviscid analysis to enable solutions to be obtained along the 
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real radial axis. With the further assumption that the mean flow properties 
are functions of the radial coordinate alone, the "parallel- flow approxima- 
tion," the linearized disturbance equations may be written, 



The velocity components in the (r, z) directions are v, w, and u, 
respectively, overbars denote the mean flow, and primes denote disturbances. 
All quantities have been nondimensionalized with respect to the jet exit 
conditions, e.g. u •= u*/uj, r ■= P ® > where stars denote a 

dimensional quantity and the subscript J denotes jet exit conditions. The 
Reynolds number and Mach number are given by, 

R *= PjUjrj/yj and M «= Uj/a^ ,• (6) 

where a is the local speed of sound. For R equations (l)-(4) approach 
the invlscid system of equations and for M = 0 they are the viscous 
incompressible equations in cylindrical coordinates. 

A disturbance is assumed to possess a wavelike form with prescribed 
azimuthal periodicity such that, for example. 
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u* (r,<t)rZ,t) - Re[u(r)exp{i(az - u)t + n4>)}] 


(7) 


where a is the complex axial wavenumber » o) Is the real radian frequency and 
n Is the azimuthal mode number. Substituting terms of the form (7) for the 
velocity and pressure fluctuations In equations (l)-(5) yields s system of 
four coupled ordinary differential equations for the four unknown functions^ 
u(r), v(r), and p(r). These equations may be written. 


-inpv + P* “ I fv" + ^ , 


( 8 ) 


^ ^ {<?” + ~ - [a^ + + !%<?} . 


(9) 


-Jflpu + pu'v + lap “ ■“ {u" + ~ - [a^ + » 

K XT 


( 10 ) 


an('l -IM^fZp + {v* + ~ + aafi} “ 0 , (11) 

where primes denote differentiation with respect to r and, 

0 ) “ au (12) 

The boundary conditions on the fluctuations which are kinematic in origin 

13 

at r = 0 were developed by Batchelor and Gill 

V, w, u, p -»■ 0 as r 
u(0) * p(0) - 0 , n 9 ^ 0 , 
v(0) * w(0) =0 , n 9 ^ 1 , 


v(0) + iw(0) = 0 , 


n = 1 . 


(13) 
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The numerical solution of equations (8) -(11) depends on a knowledge of 
the aolutions at the edges of the integration region where the mean flow 
properties are taken to be constant. These asymptotic solutions are 
developed in the next section. 

Asymptotic Solutions 

In oi'-der to determine the form of the solution in regions of constant 
mean flow properties a series solution is developed for small r. The form 
of this solution is found to indicate the general solution of the equations. 
.For n ■> 0 equation (9) is uncoupled and equations (8), (10), and (11) 
provide a system for u, v, and p. Since this case is algebraically simpler 
than that for asymmetric disturbances its solution is given more fully and 
the corresponding solution for n 0 will only be stated. 

Consider a region in which u “ u^, p « and ® where the subscript 
o indicates a constant mean flow property. Then equations (8), (10), 
and (11) reduce to, 

+ p' = i + T ■ 

+ lap “ ^ [u" + — a^u] t 0-5) 

and -iM n p + ^ ■ + iau = 0 • (16) 

o*^ r 


From the boundary conditions (13) the form of series expansions, for small r, 
for the dependent variables is found to be. 
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i 


V - V^r + 

u » U r + U«r^ + U.r^ + — 
0 z ** 

p - P^r + **' 


After considerable algebra it Is found that 




where 


_ (a^R - ) 

A*“ 2 ^ » 

(R - IM ) 


( 17 ) 


(18) 


(19) 


and - ip fi R . (20) 

■ o o 

However, equation (18) is the series expansion for the Modified Bessel 
function of the first kind, so that. 


The corresponding forms for u and v are then readily shown to be, 

(R - iM^n ) 

a . P c (22) 

o 

(R - IM^fi ) 

’ “ -’’o “* -p - n R ' ° ' 

*^0 0 

Now equations (14)-(16) are a fourth order system of equations, thus two 
linearly independent solutions which satisfy the boundary conditions should 
exist. However, the form for p contains only one arbitrary coefficient. 
This results from p only satisfying an ’’inviscid” equation. The "viscous" 
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solution for $ U zero. With p « 0 in equations (1A)-(16) they are 
identical to the Incompressible viscous equations. Thus the “viscous" 
solutions for the velocity fluctuations are identical to their incompressible 
form given by Morris^®. The general solution to the stability equations 
may now be written down. Their form is given for large r, outside the jet 
flow, where the Modified Bessel functions are replaced by Hankel functions, 
The asymptotic solutions are, 

n 0 

P - (iX*r) , (24) 

1 o 

Ot(R — IM ) y.v 

u - A. r- n -p— (lX*r) + (iXr) , (25) 

X P K O JO 

O O 

(R - 1m\) ,, . , 

« - -A IX* ^ (lX*r) - A 3 ^ (iXr) . (26) 

^00 

(27) 
CM) 


n 0 

p n R 

P = A, ° ° j ■ (iX*r) 

^ a(R - iM n) 

u *= a; (iX*r) -f A. (iXr) 

In J n ' ' 


« " *1 "n“ + ^3 f "nil S ^ 
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Tlie real parte of X* and 'X are taken to be positive bo as to satisfy the 

outer -boundary conditions, Similar forms, in terms of Modified Bessel 

functions may be written for the solutions in the potential core or very 

close to the jet axle. It should be noted that for M * 0, X** * ff and the 

viscoue incompressible solutions are recovered and for 1/R ■* 0, 

X*^ * - M^n^o , and the compressible, inviscid solutions are found, 

0 0 

Calculations 

Mean Velocity and Temperature Profiles 

The mean velocity profile considered in the calculations is representa- 
tive of jet profiles towards the end of the potential core. The same profile 

6 10 

has been used by Michalke and Morris , The velocity profile may be 
written as, 

u ■ {l + tanh [■— - r)j} , (31) 

wheroi 0 is the local momentum thickness and its variation was used by 
10 

Michalke to represent the Influence of axisyiranetry on the jet stability. 

The mean density in the jet is related to the velocity using a Crocco 
relationship of the form, 

p « T*/|l + (T* - l)u + ~ (Y “ 1) mV uil - n)} , (32) 

Which is identical to the form used by Michalke^^ except that T’*' is defined 

as the jet static temperature ratio, 

T* ^ T,/T , (33) 

u o 

where is the ambient temperature. 



9 


Nuaerlcal Procedure 

Calculations have beer perfonned for 6 * 0,16. The numerical solution 
is started at r ■ 4.0848 where u * 6.0 x lO”*^ using the linearly independent 
solutions given by equations (24)-<30), The stability equations arc 
integrated numerically using a fixed step-size Runge-Kutta integrator with 
the linear independence of the solutions being preserved by an ortho- 
normalization procedure, Ihe numerical integration is stopped at r •* 0.2448, 
where u * 0.99999 and the numerical solutions are matched with the known 
form of the solutions in the potential core. This matching procedure leads 
to a determinant whose magnitude must be minimized to obtain the eigen- 
solution. An iterative technique based on an inverse Lagranglan interpola- 
tion procedure is used to detenaine an eigenvalue. Further details of the 
numerical procedure are given by Morris^^. The inviscid calculations, 
which are shown for comparison were obtained by solving the second order 
ordinary differential equation for the pressure on the same integration 
contour. 

Results 

In spite of the large number of parameters that enter the problem in 
compressible viscous analyses only a limited number of calculations are 
presented. The trends that were previously observed in the limiting cases, 
inviscid compressible flow and viscous incompressible flow, are also found 
in the present results. All calculations have been performed for a jet 
static temperature ratio, T* = l.Q. 
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Effect Mach Number , The variation of the local growth rate, -a^G, 

for the axisynmetrlc, n ■ 0, node and the asymmetric, n ■ 1, node, are 

shown in Figures 1 and 2, respectively. Also sho?^ is the Inviscid solution 

for M ■ 0. The results are for a local Reynolds number, R0 ■ 80. in both 

cases as the Mach number increases so the maximum growth rate decreases and 

the range of amplifying frequencies decreases. The more rapid decrease in 

amplification rate with Mach number is for the n » 0 mode. This indicates 

that for higher Mach numbers the n ■ 1 mode will be more likely to occur 

2 

naturally. This was observed in the measurements of McLaughlin et al. . 

The frequency for maximum growth rate decreases with increasing Mach number 
for both modes but is lower for the n * 1 mode. The variation of the phase 
velocity, given by oj/ct^, is shown in Figures 3 and 4 for the n “ 0 and 
n •= 1 modes, respectively. At low frequencies the phase velocity for the 
n «= 0 mode is greater than the jet centerline velocity. Increasing the 
Mach number results in a decrease in the phase velocity for most of the 
amplifying frequencies. At high frequencies where the disturbances are 
damped increasing the Mach number results in a slight increase in phase 
velocity. The variation of the phase velocity with Mach number is much 
less for the n = 1 mode, the phase velocity being approximately 0.6 times 
the jet centerline velocity for all except the lowest frequencies. The 
variation of maximum growth rate and phase velocity at the frequency for 
maximum growth for the two modes is given in Table 1. In spite of the wide 
variation in phase velocity for different frequencies for the n = 0 mode 
its phase velocity at the frequency for maximum amplification is virtually 
independent of Mach number. 
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Effect of Reynolds Number . The variations with Reynolds number and 
frequency of the local growth rate, -a^6, for the n ■ 0 and n • 1 inodes 
are sho\m in Figures 5 and 6, respectively. The results are for a Mach 
nund>er of 1,0. The inviscld result at this Mach number is also shown. 

The trends are the same as those .calculated by Morris^^ for the incompressible 
case. Increasing the Reynolds number increases the growth rate at all 
frequencies. Calculations have not been performed for damped disturbances 
with 8 uper 80 P'',c phase velocities. The existence of these modes has been 
questioned by Lessen et al,^^ and even if they exist physically it appears 
that the present formulations of stability theory are inadequate for their 
correct treatment. There is little variation with Reynolds number of the 
phase velocities for either the n « 0 or n “ 1 mdde. This is shown in 
Figure 7 where the growth rates and phase velocities, relative to the 
inviscld values, are shown for both modes as a fvmction of Reynolds number 
for a fixed frequency, o)6 ■ 0.14, which is close to the maximum amplifica- 
tion frequency for the n «= 0 mode. For 0)6 = 0.14 and M = 1.0, the local 
growth rate in the inviscld limit for the n = 0 mode is - 0.0372 and 

the phase velocity is 0.7150. For the n *= 1 mode “ 0.0458 and the 

phase velocity is 0,5993, At a local Reynolds number R6 = 640 the viscous 
growth rate is 96 percent of the inviscld value for the n ~ 0 mode and 
97 percent of the inviscld value for the n == 1 mode. It is clear that the 
inviscld stability characteristics are recovered for any appropriate ’’local” 
Reynolds number greater than, say, 500.0. At lower frequencies the inviscld 
limit will be approached more slowly since it is the factor (aR) that enters 
the equations. Thus for lower frequencies a higher value of Reynolds 
number is required to come close to the inviscld solutions. 
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‘ Conclusions 

Calculations have been presented for the viscous stability character- 
istics of compressible axisynmetric jets* Asymptotic solutions to the 
compressible stability equations have been obtained and have been used as 
boundary conditions for the numerical integration. Increasing the jet Mach 
number is found to be stabilizing as is decreasing the Reynolds number. 

The helical, n - 1, mode, is more unstable than the axisymmetric, n « 0, 
mode at high Mach numbers. The viscous calculations have been found to 
approach the inviscid limit for local Reynolds numbers greater than 500. 

A number of simplifying assijmptions have been made in the development 
of the compressible stability equations. Extension of the method presented 
here to a system of equations which includes such effects as temperature 
dependent molecular properties appears reasonable. However, if the Reynolds 
number is based on an eddy viscosity the present system of equations appears 
adequate. 

Acknowl edgmen t 

This ,:ork was supporti^'d by the NASA Langley Research Center under 
Grant NSG 1580. The technical monitor is Dr. T. B. Gatski, 



13 


References 

^ Moore, C. J., "The Role of Shear Layer Instability Waves in Jet 
Exhaust Moise." Journal of Fluid Mechanics , Vol, 80, Pt. 2, 1977, 
pp. 321-367. 

^ McLaughlin, D. K., Morrison, G. L., and Troutt, R. R. , "Experiments 
on the Instability Waves in a Supersonic Jet and Their Acoustic Radiation," 
Journal of Fluid Mechanics , Vol. 69, Pt. 1, 1975, pp. 73-95. 

■3 Tam, C. K. W. and Morris, P. J., "The Radiation of Sound by the 
Instability Waves of a Compressible Plane Turbulent Shear Layer," Journal 
of Fluid Mechanics . Vol. 98, Pt. 2, 1980, pp. 349-381. 

^ Bechert, D. and Pfizenmaier, E., "On the Amplification of Broadband 
Jet Noise by a Pure Tont Excitation," AIAA Paper 76-489, July 1976, 

^ Lessen, M., Fox, J. A., and Zien, H. M., "The Instability of 
Inviscid Jets and Wakes in Compressible Fluid," Journal of Fluid Mechanics , 
Vol. 21, Pt. 1, 1965, pp. 129-143. 

^ Michalke, A., "Instabilitat eines kompressiblen runden Freistrahls 
unter Berucksichtigung des Einflusses der Strahlgrenzschichtdicke, " 

Z. Flugviss . Vol. 19, 1971, pp. 319-328. 

^ Burridge, D. M,, "The Stability of Round Jets," Ph.D. Thesis, 

Bristol University, 1968. 

8 

Mollendorf, J. C. and Gebhart, B., "An Experimental and Numerical 
Study of the Viscous Stability of a Round Laminar Vertical Jet With and 
Witho'ut Thermal Buoyancy for Symmetric and Asymmetric Disturbances," 
Journal of Fluid Mechanics, Vol. 61, Pt. 2, 1973, pp. 367-399. 



Leesen, M. and Singh, P. "Tlie StahlXity of Axisyimnetrlc Free 
Shear Liiyers," Journal of Fluid Mechanics , Vol, 60, Ft. 3, 1973, pp. 433-457. 

HorrlB, P. J., "The Spatial Viscous Instability of Axlsynnoetric 
Jets," Journal of Fluid Mechanics , Vol. 77, Pt. 3, 1976, pp, 511-529. 

Lees, L. and Lin, C. C., "Investigation of the Stability of the 
Laminar Boundary Layer in a Compressible Fluid," NACA Tech. Note No. 1115, 
1946. 

12 

Dunn, D. W. and Lin, C, C., "On the Stability of the Laminar 
Boundary Layer in a Compressible Fluid," Journal of the Aeronautical 
Sciences . Vol. 22, 1955, pp. 455-477. 

13 

Batchelor, G, K, and Gill, A. E., "Analysis of the Stability of 
Axisymirtetric Jets," Journal of Fluid Mechanics , Vol, 14, Pt. 4, 1962, 
pp. 529-551. 

14 

Lessen, M., Fox, J, A., and Zien, H. M. , "Stability of the Laminar 
Mixing of Two Parallel Streams with Respect to Supersonic Disturbances," 
Journal of Fluid Mechanics , Vol. ,25, Pt. 4, 1966, pp. 737-742. 



15 


Table 1. Variation of naxinun amplification rate and phase velocity at the 
maximum amplifying frequency with Mach number, R6 - 80. 



n “ 0 

n « 1 

M 

max. 

amplifying 

frequency, 

0)6 


o)/a 

r 

max. 

amplifying 

frequency, 

0)0 


o)/a 

r 

0 

0.1698 

O.O490 

0.7226 

0.1278 

0.0515 

0.6051 

0.4 

0.1638 

0.0454 

0.7230 

0.1236 

0.0497 

0.6021 

0.8 

0.1474 

0.0351 

0.7224 

i 

0.1110 

0.0447 

0.5930 

1.0 

0.1348 

0.0278 

0.7228 

0.1026 

0.0413 

0.5862 

1.2 

0.1201 

0.0188 

0.7225 

0.0928 

0.0375 

0.5780 
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Fig. 1 


•Fig. 2 


Fig. 3 


Fig. k 


Fig. 5 


Fig. 6 


Fig. 7 
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Figure Captions 

Effect of Mach number on the variation of amplification factor 
with frequency, n » 0, R6 « 80,0. M ■ 0, — — — — ; 

M « 0,4,— ' — ; M ■ 0.8,-- ; M « 1.0, 

M ■ 1.2,— ——— . Inviscid solution, M * 0, , 


Effect of Mach number on the variation of amplification factor 
with frequency, n * 1, R6 « 80.0. For legend see Fig. 1, 

Effect of Mach number on the variation of phase velocity with 
frequency, n ** 0, R6 - 80.0. For legend see Fig. 1 

Effect of Mach number on the variation of phase velocity with 
frequency, n = 1, R0 = 80.0. For legend see Fig. 1. 


Effect of Reynolds number on the variation of amplification 
factor with frequency. n*=0, M = 1.0. R0 = 20.0,-' — 

R6 = 40.0, . R0 = 80.0, — — ; R0 = 160.0, 

— — — ; R0 = 320.0, — — — — — , Inviscid solution, 

M = 1.0, . 

Effect of Reynolds number on the variation of amplification 
factor with frequency, n = 1, M = 1.0. For legend see Fig. 5. 
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